From quantum theory to classical dynamics under spontaneous 

wave-packet reduction 



C. F. Huang 

1 Department of Physics, National Taiwan University, Taipei, Taiwan, R. 0. C. 
2 National Measurement Laboratory, Center for Measurement Standards, Industrial Technology 
Research Institute, Hsinchu, Taiwan 300, R. 0. C. 
(February 1, 2008) 



Abstract 



A quantum master equation of the Lindblad form is obtained in this pa- 
per by considering the spontaneous wave-packet reduction. Different classical 
equations can be derived after exactly mapping such a quantum master equa- 
tion to a continuous time random walk (CTRW). Although this CTRW is a 
quasiclassical walk, the effects due to the quantum interference can still be 
important in such a walk. Macroscopically, we shall consider the uncertainty 
of the potential and determine the effective transition probability by a family 
of Schrodinger equations (or operators). 



I. INTRODUCTION 



While classically a particle has the definite position and momentum at any time, in quan- 
tum mechanics it is described by a complex wavefunction and its position and momentum 
cannot be determined simultaneously from the uncertainty principle. [1,2] A wavefunction, 
in fact, can be a superposition of microscopic wave packets and become very extended with 
respect to either the momentum or position. To determine the momentum (position), in 
quantum mechanics it is assumed that the wavefunction collapses in momentum (position) 
space under the momentum (position) measurements. [1,3] With suitable assumptions, the 
quantum master equation of the Lindblad form [4,5] can be used to explain the appearance 
of the classical world. [6] In addition to the quantum Liouville term, in such an equation 
there is a relaxation term which can induce the collapse of states or reduction of wave packets 
[1,3,6-8]. The relaxation term can be due to the quantum measurements or the interactions 
with the reservoir. [6,7,9-11] Although no state has the definite momentum and position 
simultaneously, a wavefunction can collapse to a wave packet to have the momentum and 
position within a small uncertainty under the phase-space measurement [9,12]. A classical 
particle, in fact, is described by a microscopic wave packet of which the uncertainties in 
both the momentum and space are negligible macroscopically. [1] Let |0) be a normalized 
microscopic wave packet located at the position x = with the momentum p = 0. Then in 
the set 



each ket is a wave packet located at x = x with p = p . Here T xo is the translation operator 
to shift each ket by the displacement x , and T po is the operator to shift each ket by the 
momentum po- 

In this paper, a particle moving in R n space under the Hamiltonian 




(1) 




(2) 
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is considered, where n is an integer, m is the mass of the particle, and V(x) is the potential. 
Let p{t) be the density matrix at time t and satisfies the normalized condition 

trpit) = 1. (3) 

(All the quantities and equations are discussed in the Schrodinger picture [1] in this paper.) 
Considering a lifetime r for the particle to become wave packets in S<f, spontaneously, in 
section II the equation 

= jrWt),m + ^pi^p)(0x,p|p(oi0x,p)(0 x ,p|. (4) 

is derived under suitable assumptions, where H is the Plank constant divided by 2n. Eq. 
(4), in fact, is a quantum master equation of the Lindblad form and can be exactly mapped 
to a continuous time random walk (CTRW) [13,14] in which the probability density 

r(x,p,i) = -(0x,p|p(Ol0x, P )- (5) 

T 

In such a CTRW, r(x, p, t) is governed by 

r(x, p, t) = £ dt' J d ^ n y ^(x, p; x', p'; t - t')r(x', p', t') + s(x, p, t). (6) 
with the transition probability density (memory kernel) 

^(x,p;x',p';f-0 = ^e-^'^K^ple-t^'^l^^)! 2 (7) 
and the source term 

S (x, P ,t) = i e -^)/M0x, P |e-^ t - t <^p(to)e^ t - t <^|0 x , P ). (8) 

Here to is the initial time, and in Eq. (6) / d ^ 1 ^ l y presents the sum over all possible states 
in the CTRW with (27r^) _ ™ as the density of states in phase space. Since r(x, p, t) is a real 
function defined in the (classical) phase space at any time t, such a CTRW can be taken 
as a quasiclassical walk. As shown in section III, however, the effects due to the quantum 
interference can be important even in the macroscopic scale. 
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In the quantum master equation of the Lindblad form, the irreversibility is due to the 
relaxation term. On the other hand, E. T. Jaynes [15] obtained a theory for irreversible pro- 
cesses by considering a family of Schrodinger equations (or operators) under the assumption 
that the potential may not be known certainly. The random Schrodinger equations (or oper- 
ators) [16-18], in which the potentials are random fields, have been used to study phenomena 
due to the quantum interference such as the Anderson localization [18,19]. Macroscopically, 
as shown in section III, we shall consider the uncertainty on the position of the potential 
and determine the effective transition rate by a family of the Schrodinger equations (or op- 
erators). With suitable assumptions, the CTRW can correspond to the classical Liouville 
equation [20] or classical master equation [10,20]. 

Different quantum approcahes have been developed to derive the classical equations. 
[6,15,16,20-24] It is shown in section IV that Eq. (4) can be related to the momentum- 
position collapse model discussed in Refs. [23] and [25]. The density ^(x, p; x', p'; t — t') 
given by Eq. (7) is the product of a temporal factor and the well-known quantum transition 
probability for being |</> x ',p') at time t! to become |0 X)P ) at time t. Therefore, the CTRW 
provides an intuitive way to relate the quantum theory to the classical motions under the 
wave-packet reduction. The transition probability density ^(x, p; x', p'; t — £'), in fact, is 
proportional to the jump rate in the corresponding quantum trajectory formulation. [26,27] 
The conclusion is made in section VI. 

II. QUANTUM MEASUREMENTS, REDUCTION OF WAVE PACKETS, AND 

CONTINUOUS TIME RANDOM WALK 

Under the (quantum) Liouville flow with the Hamiltonian H, the time-dependent density 
matrix satisfies [1,11] 

pit + At) = e~* AtH p(t)ei AtH . (9) 

But if we performed a measurement of the nondegenerate observable A at t + At, the particle 
must collapse suddenly to an eigenket of A at t + At. [1] Let 
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{\a)\ A\a) = a\a), (a\a) = 1}. (10) 

be the set of normalized eigenkets of A. (Note that the eigenkets of momentum or position 
operators cannot be normalized in L 2 . [1]) Since the probability to collapse to \a) is 

c a (t,At) = (a\e^ AtH p(t)e^ AtH \a), (11) 

at t + At we shall replace Eq. (9) by 

p(t + At)=J2c a (t,At)\a}(a\ (12) 

a 

under such a measurement. [1] 

Now assume that the particle has the lifetime r to become to an eigenfunction of A 
spontaneously rather than collapsing suddenly. Then during small At the probabilities to 
collapse and not to collapse are At/r and 1 — At/r, respectively, and at t + At we shall set 

pit + At) = (1 - At/r)e-^ AtH p(t)e^ AtH (13) 
+ (At/r)J2ca(t,At)\a)(a\ + o(At 2 ). 

a 

At the right hand side of the above equation, the first term presents the factor due to the 
Liouville flow and the second term presents the factor due to the collapse. As At — > 0, with 
Eq. (11) we can reduce Eq. (13) as 

|p(t) = t m ,H] - e^L + I ^ \ a )(a\p(t)\a)(a\. (14) 

It is easy to check that trp(t) is kept under the above equation. With some calculations, we 
can see that the last two terms in Eq. (14) induce the decay of phases (off-diagonal terms) 
with respect to the eigenkets of A without affecting the diagonal terms. 

Although a wavefunction cannot collapse to a state with the definite momentum and 
position from the uncertainty principle, it can collapse to a ket in S<f> defined in Eq. (1). To 
replace eigenstates of A by the states in S^, the third term of Eq. (14) should be replaced by 
°I fSpl0x, P )(0x,p|p(i)|0x,p)(0x,p| with c as a constant. (We shall replace £ a by / 
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because the momentums and positions are continuous parameters.) To preserve trp(t), we 
shall set c = 1/r and Eq. (4) is obtained. When the spatial width of \<f>) tends to zero, each 
ket in becomes an eigenket of the position operators and hence Eq. (4) describes the 
collapse in position. On the other hand, Eq. (4) describes the collapse in the momentum 
space as the width of \<f>) in momentum shrinks to zero. 

The density matrix p(t) governed by Eq. (4), in fact, satisfies 

p( t ) = e -(t-to)/T e -i(t-t )H p ^ e ±(t-t )H + (15) 

£ Mew J g^(x, p, Oe-^^p) 

with r(x, p, t) defined in Eq. (5). To obtain the above equation intuitively, note that the 
first term at the right hand side of Eq. (15) corresponds to the case that the particle 
moves without collapsing. The matrix e~^^~ to ^ H p(t )e^^~ to ^ H corresponds to the quantum 
Liouville flow and the factor e~^~ to ^ T is the probability for the particle not to collapse. 
On the other hand, if the particle collapses and the corresponding density matrix becomes 
|0x,p)(0x, P | at time t', at time t > t' the density matrix is e~^ i_ *'^|0 XiP ) (0 x , p |es(*~*') H under 
the quantum Liouville flow. The probability to become |0 X)P )(0 X , P | at t' is proportional to 
r(x, p, t') = ^(</>x, P |p(Ol0x,p) and the collapse can still occur when t > t', so we can expect 
the last term in Eq. (15). Because Eq. (4) is a first-order differential equation with respect 
to t, to prove Eq. (15) we just need to take the time derivative on the above equation and 
check the initial condition. Then inserting the above equation into the right hand side of Eq. 
(5), we can obtain Eq. (6) with the transition probability density and source term defined 
in Eqs. (7) and (8). 

As shown in Appendix A, actually Eq. (4) can be taken as a particular case of the 
following equation 

J^(0 = ^W,w(01-^E^^^ (is) 

Here Ti(t) is a time- dependent Hamiltonian, all states in the last two terms are from a 
set S composed of normalized kets, the coefficients ww are nonnegative real numbers, and 
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{A, B} = AB + BA for any two operators A and 5. The above equation is just the master 
equation of the Lindblad form 

= ^P(t),n(t)] - l^{ P (t),viv a } + ^v a p(t)vi (17) 

1 II 

with a — > (/, /') and the operator V Q — > ft is shown in Ref. [10] that Eq. (16) can 

be mapped to an extended random walk, and Eq. (15) in this paper can be obtained from 
Eq. (13) in Ref. [10], in fact. (But note that in Ref. [10] the Hamiltonian is time-dependent 
and the transition probability density is not necessary of the form l,t — t').) To prove 
that Eq. (6) defines a CTRW, we need to show that [10,13,14] 

r°° , r d n xd n p , . . . .. 

L */(^^.p;<p';*-*') = i. (is) 

r°° , r d n xd n p . , 

i dt J^w s{ ^ =1 - (19) 

The proof is given in Appendix A. 

III. FROM QUASICLASSICAL CTRW TO CLASSICAL EQUATIONS OF 

MOTION 

By considering a lifetime r for the particle to become a microscopic wave packet spon- 
taneously, in the last section Eq. (4) is derived and is mapped to a CTRW no matter how 
long the lifetime r is. In this section, it will be shown that macroscopically such a CTRW 
can be reduced to the classical Liouville equation, classical master equation, or a random 
walk exhibiting effects due to the quantum interference. (It should be emphasized that some 
assumptions used in this section do not hold in the conventional quantum scattering theory. 
[1,2]) 

Because the source term given by Eq. (8) contains a decay factor e~^~*°^ T , such a 
source term is unimportant when t » t and the statistical properties of the CTRW are 
determined by the transition probability density in Eq. (7). Define 
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G(x,p;x',p';0 = | (0 x , p |e^^|0 x ^) | 2 (20) 

as the expectation value of the time-dependent ket |0 X ', P '(^)) = e~^ tH \4>x',p') i n ket |0 X , P ) at 
t = £. We can rewrite the transition probability density as 

V(x, p; x', p'; = V^G(x, p; x', p'; (21) 

r 

with £ = t — t'. The time-dependent ket |0 x ', p '(£)), in fact, is governed by the Schrodinger 
equation 

iK§j\fa,At)) = H\<f> x >,At))- ( 22 ) 

To determine the behavior of the CTRW, we only need to discuss G(x, p; x', p'; £) or 
|0x', P '(O) when £ is comparable with r because of the factor \e~ s -l T at the right hand 
side of Eq. (21). 

Although the CTRW defined by Eqs. (5)- (8) is a walk over the phase space and can 
be taken as a quasiclassical random walk, the effects due to the quantum interference may 
still be important in such a CTRW. To see this, let V(~x) in Eq. (2) be a random potential. 
While classically a particle with nonzero velocity can move to infinity if V(~x) is very weak, 
the quantum interference can make |0 x ' iP '(t)) governed by Eq. (22) become immobile after 
passing long enough time. [19] In the one or two-dimensional cases, in fact, insulating behav- 
iors are expected no matter how weak the random potential is. [28,29] Since G(x, p; x', p'; £) 
is the expectation value of |0 x ' lP '(£)) i n |0x, P ) and the behavior of the CTRW is determined 
by G(x, p; x', p'; £) when £ is comparable with r, we can expect that the effects due to the 
quantum interference exist if r is very long. It should be noted that in Eq. (15) r(x, p, t') 
is followed by the matrix e - ^* - *'^ |0 X)P ) (0 x , p |e^*"*') H , which is an immobile state rather 
than a state with the velocity p/m when t — t' is large. Such an immobile state is the su- 
perposition of the microscopic states in S^, and can be a wave packet of macroscopic scale. 
On the other hand, consider the case that in Eq. (2) the change of V(x) is insignificant 
in the quantum scale and m is very large. In such a case, it is known that Eq. (22) can 
correspond to the classical Liouville equation for a very long time. Hence we can expect 
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that with suitable r, the CTRW corresponds to the Liouville equation in such a scale. The 
details are given in Appendix B. 

It is shown in Appendix C that under suitable assumptions, macroscopically the effective 
transition probability density of the CTRW is 

*(x, p; x', p'; t-t') = — L- / d n y f d n k^(x + y, p + k; x' + y, p' + k; t - t') (23) 

il^h ^|y|<Ai J|k|<A 2 

Here Ai and A 2 are real parameters so that two points (xi,pi) and (x 2 ,p2) in the phase 
space are too close to be distinguished macroscopically if |xi — x 2 | < Ai and |pi — p 2 | < A 2 , 
and fli and f2 2 present the volumes j| y |<Ai d n y and J| p |<a 2 d n p- That is, the macroscopic 
jump rate is the average of the microscopic transition probability on a small volume which 
is taken as a point macroscopically. Define H y = ^ + V(x — y) = T^HT y with y G R n . 
The density \l/(x, p; x', p'; t — t'), in fact, is determined by the set {H y , |y| < Ax} and hence 
can be related to the irreversible theory suggested by E. T. Jaynes in Ref. [15]. To see this, 
note that Eq. (23) can be rewritten as 

e -(t-t')/r . 

(x, p; x', p'; t-f) = / (0x, P+ kkx',p' + k(t - f) |0 XiP+k )rf n /c (24) 

if we define 

^x, P (0 = J <TyF(y)U<y\Z)aWuW(Z). (25) 

with a^ p = |0 x ,p)(0x,p|- Here the real function F(y) equals il^ 1 for |y| < A 1 and equals 
for |y| > Ai, and f/ (y) (0 = e"*^ with f real parameter. Eq. (25) can correspond to 
Eq. (12.5) in Ref. [15] after taking F(y)d n y as the probability to move the matrix a^ p by 
the unitary operator U^(^). To prove Eq. (24), note that from the definition of 

|0x +y , P )=^x +y T p |0)=T y |0 XiP ) (26) 

and from the Taylor's expansion 

e-iW* = £ ^(~) j (t - t'y(T y HT y y (27) 
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Since each H y is obtained from H by performing the position translation on the potential, 
here the randomness of the time evolution comes from the uncertainty of the potential in 
the macroscopic scale. 

As mentioned above, with suitable assumptions on r, m, and V(x), the transition prob- 
ability density ip and the statistical properties of the CTRW can be determined by the 
classical Liouville equation. If V(x) is of the macroscopic scale so that V(x.) ~ V(x + y) 
when |y| < A 1; it is easy to see that macroscopically the derived CTRW still corresponds to 
the classical Liouville equation. But if V(x) is random in the scale of A 1; in the two or three 
dimensional cases the CTRW may correspond to the classical master equation in a large 
scale since classically the (microscopic) Liouville equation can be reduced to the classical 
master equation. [20] On the other hand, in the one or two dimensional cases insulating 
behaviors due to the quantum interference are still expected macroscopically when V(x) is 
random and r is long enough for the quantum interference. Therefore, macroscopically the 
CTRW defined by Eqs. (5)-(8) may correspond to the Liouville equation, master equation, 
or a random walk exhibiting the effects due to the quantum interference. 



Different quantum approaches have been developed to obtain classical equations. 
[6,15,16,20-24] In this paper, Eq. (4) is derived by considering a relaxation time for wave- 
functions to collapse to microscopic wave packets, which are concentrated in both momentum 
and position. In Refs. [23] and [25], a spontaneous momentum-position localization model 
is discussed based on the equation of the following form 



Here a and (3 are two positive real parameters, r is the relaxation time, and c is a parameter 
which can be obtained by considering the conservation of trp(t). (Thus c is determined by 



IV. DISCUSSIONS 



d_ 
di 




(28) 
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a, 0, and r.) Let Sj£ = a(p - k) 2 + /3(x - y) 2 . We can rewrite Eq. (28) as 

Define /i"'^ = ap 2 + /3(x — y) 2 as the Hamiltonian of a simple harmonic motion. Then Sy'£ 
can be obtained by shifting hy' 13 with the momentum translation operator. [23] Thus we 
can obtain the eigenkets and eigenvalues of S^'f and exp(—Sy'yj exactly. Let |-Ey '£/ be the 
ket satisfying S*{\Kt) = KilKi)- E Q- ( 4 )> in fact > can be obtained from the above 
equation by considering the limit 

a — > oo and /3/a = const. (30) 
In such a limit, each \Eyfy is unchanged while Ey'£ — * oo. And it is easy to see that 

exp(~S%) = E e- E "AKi)( E yi\ ( 31 ) 

where £y^ is the lowest eigenvalue of Sy 'f. Each ket |£y'k)> in fact, is just a Gaussian wave 
packet, and Eq. (4) is obtained after setting \<f>) as a wave packet of Gaussian type. 

The quantum trajectory formulation is a powerful numerical approach to solve Eq. (17), 
the master equation of the Lindblad form. [26,27] In such a formulation, stochastic quan- 
tum jumps interrupt the continuous evolution determined by the nonhermitian effective 
Hamiltonian, 

% = ^-S^ (32) 

z a 

The continuous evolution is responsible for the first two terms at the right hand side of Eq. 
(17) while the last term in such an equation is interpreted as the origin of quantum jumps. 
Each V a is the "collapse" operator to induce the collapse of an arbitrary ket \a) to V a \a) in 
a quantum jump. For Eq.(4), the effective Hamiltonian is reduced as 

H eff = H - U/2. (33) 
11 



In addition, we can set a = (x, p) and take each 

^|0x,p)(0x,p| (34) 

as a collapse operator. In a trajectory, a wavefunction becomes a wave packet at a quantum 
jump when each collapse operator is given by the above equation. After the first jump, it 
is easy to see that the following (conditional) jump rate is proportional to the density ip 
given in Eq. (7) with H e ff defined in Eq. (33). Therefore, we just need to consider the 
jump rate proportional to the transition probability density for the CTRW to determine 
statistical properties. 

We can see that r(x, p, t) for the discussed CTRW is positive definite from its definition 
while the Wigner distribution [20], which is taken as the classical correspondence of p(t), 
may contain the negative part. The density r(x, p,t), in fact, only records the probability 
for a wavefunction to become to a wave packet at time t. Therefore, its physical meansing 
is different from that of Wigner distribution. 

V. CONCLUSION 

In this paper a quantum master equation of the Lindblad form is obtained by considering 
the spontaneous wave-packet reduction in phase space. Such an equation can be mapped 
to a quasiclassical continuous time random walk, from which the classical master equation, 
classical Liouville equation, or a walk exhibiting effects due to the quantum interference are 
obtained in the macroscopic scale. The macroscopic transition rate, in fact, is determined 
by a family of Schrodinger operators. 
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APPENDIX A 



In Eq. (16), the transitions are over a set of normalized kets and it is not necessary to 
ask that such a set is composed of orthogonal kets. [10] Setting 



w U ' = 6 U '/t, (35) 



Eq. (16) can be reduced to 



j t P{t) = {W)^)] - Y T {p{t)^\l){l\} + \ £ \l)(l\p(t)\l)(l\. (36) 

In the case that all transitions are over the set S^, in the above equation each / should 
correspond to a wave packet XiP and we shall replace J2i by / d ( 2 wh)™ ■ After some calculations, 
we can obtain 

/w luw = ' (37) 

and then derive Eq. (4) from Eq. (36) by setting H(t) = H, where I is the identity operator. 
Therefore, Eq. (4) is just a particular case of Eq. (16). 

To prove Eqs. (18) and (19), note that for any operator O, 

f fe"^/ ^tr(Oe^-')-|0 x , p )(0 x , p |e-(-')-) (38) 



= r * e -M)/rw 0e iM)fl e -i('-^)) (from Eq. (37)) 

Jt' T 

= trO. (since e"s (t ~*' )H is unitary) 

From the above equation we can prove Eq. (18) by taking O as |0 x ',p')(0x',p'|, an d prove 
Eq. (19) by setting O and t' as p(t ) and to, respectively. 
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APPENDIX B 



Consider the case that the change of V(x) is insignificant in the quantum scale and m is 
very large. For the particle corresponding to the wave packet |0 X ', P ') initially, it is known that 
the quantum effect can be ignored for a long time. Hence the function (j(x, p; x', p'; £), which 
is the probability to find this particle around (x, p) at time £ under the Liouville flow, con- 
centrates around the point (x, p) = (X(x', p', £), P(x', p', £)). Here (X(x', p', £), P(x', p', £)) 
is the location of the particle which originally is at (x', p') after passing time £ under the 
classical Liouville flow. 

Approximating G(x, p; x', p'; £) as 5(x — X(x', p', £))<5(p — P(x', p', £)) in Eq. (21), from 
Eq. (6) it is not hard to see that (0 x ,p|p(£)|0x,p)(= r(x, p, t)r) satisfies 

|(0 x ,p|p(t)l0x, P ) + i ^^-(0x, P |p(t)l0x,p) - 1 ^^:(0x,pIpWI0x,p> - o 

3 — 1 ^ 3 — 1 ^ ^ 

as i >> t - Taking (0 x , P |p(£)|0x,p) as the distribution function, therefore, the CTRW defined 
by Eqs. (5)- (8) can be related to the classical Liouville equation. 

APPENDIX C 

Let 

1 r t+T r r 
i?(x, p, t) = — — / dt' / d n y / d n k r(x + y, p + k, t') (39) 

\ll\l2l Jt J|y|<Ai J|k|<A 2 

be the macroscopic probability density of the CTRW. Here T is a real parameter so that for 
any two times t\ and ti satisfying \t\ —ti\ < T, they are indistinguishable macroscopically. 
Assume that microscopically r(x, p, t) is homogeneous in the time scale of T so that 

1 r t+T 

- J t dt' r(x + y, p + k, t') ~ i?(x, p, t) (40) 

when |y| < Ai and |k| < A 2 . To show that \l/(x, p; x', p'; t — t') in Eq. (23) can be taken as 
the macroscopic transition rate of the CTRW, we need to prove that in a large scale 

i?(x, p, t) ~ J £ df*(x, p; x', p'; t - t')^(x', p', + 5(x, p, f) (41) 
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under Eq. (40), where S(x, p,t) is the macroscopic source term. To prove that the above 
equation holds in a large scale, we just need to insert Eq. (6) into Eq. (39) and note Eq. 
(40). 
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